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Abstract: This paper includes a proof of well-posedness of an initial-boundary value problem involving 
a system of non-local parabolic partial differential equation(PDE) which naturally arises in the study of 
derivative pricing in a generalized market model which is known as a semi-Markov modulated geometric 
Brownian motion(GBM) model. We study the well-posedness of the problem via a Volterra integral equation 
of second kind. A probabilistic approach, in particular the method of conditioning on stopping times is used 
for showing the uniqueness. 
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1 Introduction 

In the Black-Scoles-Merton model of option pricing the dynamics of stock price {St}t>o is assumed to be 
given by a geometric Brownian motion (GBM), i.e., 

dSt = ^iStdt aStdWt, Sq > 0 

where /r and cr are positive constants and {Wt}t>o is a standard Wiener process. Existence of a deterministic 
risk free asset with constant growth rate r is also assumed. Under these model assumptions, the European 
call option price function satisfies a parabolic partial differential equation 

d d 15^ 

—V?(f, s) + s) + s)=r s). 

This is known as Black-Scholes-Merton partial differential equation (B-S-M PDE), whose solution can be 
obtained explicitly [3]. However, empirical studies of financial assets suggest that the assumption of constant 
or deterministic cr, the volatility coefficient, is rather unrealistic. To overcome this drawback, there has been 
an increasing interest in market models where cr is taken as a Markov process. In stochastic volatility 
models, the square of cr is taken as an Ito diffusion |14) . whereas in regime switching models cr is driven by a 
finite-state continuous-time Markov chain m. i, 0 , 0 , m. m. na, m and m- The resulting market 
models for both of theses types are incomplete. Option pricing in such markets has been studied by several 
authors using different approaches. Using locally risk minimizing approach, the European call option price 
in a regime switching market is shown to satisfy a generalized B-S-M PDE 0 , 0 , 0 - 

Locally risk minimizing option pricing in a regime switching model with semi-Markov regimes is studied in 
0. There the market parameters, /i, cr and r are driven by a finite-state semi-Markov process {Xt}t>o- 
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Sojourn or holding times in each state in a finite-state continuous-time Markov chain are distributed as 
exponential random variables, whereas that could be any positive random variable for semi-Markov case. 
Thus the class of semi-Markov processes subsumes the class of Markov chains. There are some statistical 
results in the literature (see [15] and the references therein) which emphasize the advantage of use of semi- 
Markov switching models over simple homogeneous Markov switching models. For example it is mainly 
useful to deal with the impact of a changing environment (i.e. the business cycle), which exhibits duration 
dependence. Suppose {Xt} is a semi-Markov process and the stock price {S'*} is given by 

dSt = ^i{Xt)Stdt +(j{Xt)StdWt, So > 0. 

Then it is shown in that the call option price function satisfies a non local system of parabolic PDF 

d d 9 19^ 

—(^(t, s, i, y) + s, i, y) + r{i)s—if{t, s, i, y) + s, i, y) 

+ 1 = r{i) ip{t,s,i,y), 

1 - ^(2/ I ^ 

defined on 

V := {(t, s, i, y) e (0, T) x R+ x T x (0, T) \ y G (0, t)}, (1) 

and with conditions 

\imip{t,s,i,y) = 0 Vte[0,r] 

siO 

^{T,s,i,y) = {s-K)+; se(0,oo); 0<y<T; i = l,2,...,fc (2) 

where k is the number of all possible states of Xt, r{Xt) is the spot rate of interest at time t, (pij) are the 
transition probabilities to state j form i, F{- \ i) is the conditional distribution function of holding time, 
f{y I i) = -^jjFly \ i), K is the strike price and T is the maturity time. 

In this paper we replace the semi-Markov process {Xt} by a more general class of age-dependent processes 
which is a much wider class of switching than that appear in [9]. An age-dependent process {Xt}t>o is 
specified by its instantaneous transition rate which is a collection of measurable functions Xtj : [0, oo) —>• 
(0, oo) where i j G X := {1,2,... ,k}. Indeed {Xt}t>o is defined as the strong solution of the following 
system of stochastic integral equations 

Xt = Xo+ [ [ /i(A„_,y„_,z)p(dM,dz) (3) 

d(0,t] JR 

Yt = t- / g{Xu-,Yu-,z)p{du,dz) 

d(0,t] JR 

where p{dt,dz) is the Poisson random measure with intensity dtdz, independent of Xq and 
Hhy,z)-= O'-0lAy(y)(2:), 9{i,y,z)-- 

j £ A' i £ A’ 

where X.ij{y) are the consecutive (with respect to the lexicographic ordering on X x X) left closed and right 
open intervals of the real line, each having length Xij{y). We refer [10] for more details about this kind of 
pure jump processes. Under some smoothness and tail assumptions on Xij{y) and independence of W and 
p, we obtain the following equation of locally risk minimizing price of call option 

9 9 9 1 9^ 

—p(t, s, i, y) + s, i, y) F r{i)s—if{t, s, i, y) + s, i, y) 

+ '^^ijiy){‘P{t,s,j,0) - (p{t,s,i,y)) = r{i) (p{t,s,i,y), (4) 

defined on TX as in o and with conditions as in (ED- Some of the special cases of this equation appear in 
i, 0, US], [9], |18j . and |6] for pricing a European contingent claim under certain regime switching market 
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assumptions. Owing to the simplicity of the special cases, generally authors refer to some standard results in 
the theory of parabolic PDE for existence and uniqueness issues. But in its general form which arises in this 
paper, no such ready reference is available. In this connection we would like to emphasize that o is a non¬ 
local degenerate parabolic PDE on a non-rectangular domain. Therefore, we produce a self contained proof 
using Banach hxed point theorem. These we accomplish in two steps. First we consider a Volterra integral 
equation(IE) of second kind and establish existence and uniqueness result of that. Then we show in a couple 
of propositions, that the PDE and the integral equation(IE) problems are “equivalent”. Thus we obtain the 
existence and uniqueness of the PDE in Theorem 14. 2 1 Some further properties, viz the positivity and growth 
property are also obtained. It is also shown here that the partial derivative of the solution constitutes the 
optimal hedging strategy of the corresponding claim. We further show that the partial derivative of ip, can 
be written as an integration involving p which enables one to develop a robust numerical scheme to compute 
the Greeks. This study paves the way for addressing many other interesting problems involving this new set 
of PDFs. 

The rest of this paper is arranged in the following manner. The market model assumption and the pricing 
approach is described in Section 2. In Sections 3 and 4 the well-posedness of the Cauchy problem Qis 
established. In Section 5, the well-posedness result is used to solve the pricing and hedging problem. We 
conclude our paper in Section 6 with a few remarks. 


2 Market model 

We consider the maps Xij : [0,oo) —>■ (0,oo) where i ^ j G X := {l,2,...,fc} and define \ii{y) := 
“ Sj/zjgAf € X and y G [0,oo). We further assume the following. 

(Al) (i) For each i and j, A^- is a differentiable function, 

(ii)limj,^oo Ai(y) = oo, where Ki{y) := /J' \Xii{y)\dy. 

Define F{jj | i) := I - exp(-A,(y)), f(y \ i) := ^F{y \ i) and for each i ^ j, p^j{y) := with pu{y) = 0 

for all i and y. Set ptj := pij{y)dF{y \ i). 

(Al) (iii) The matrix {pij) is irreducible. 

Let (D, F, P) be the underlying complete probability space where a Poisson random measure p and a standard 
Wiener process W are defined and they are independent. Let {Xt,Yt) be the strong solution of system of 
equations ([3]) where the maps Xij satisfy (Al) (i)-(iii). From ([3]) it is apparent that Xt is a right continuous 
(since the integrations are over (0,t]) jump process having left limits and taking values in X. Let T„ denote 
the time of nth jump of A*, whereas Tq := 0 and r„ := Tn — Tn-i- For a fixed t, let n{t) := max{n : T„ < t}. 
Thus Tn(t) < t < and Yt = t — T^^t). It is shown in [TU] that F(y \ i) is the cumulative distribution 

function of holding time and Pij(jj) is the conditional probability that X transits to j given the fact that it 
is at i for a duration of y. From (Al)(ii), limj^_>oo F{y I *) “ 1- Hence, sojourn times are finite almost surely. 
We also note that for i ^ j, Xij{y) = hold. Thus the instantaneous transition rate is given by 

the maps Xij. 

Let {Hjjoo be the price of money market account at time t where, spot interest rate is = r{Xt) and 
Bq = 1. Thus we have Bt = e-^oLet be the price process of the stock, which is governed by 

a semi-Markov modulated GBM i.e., 

dSt = y{Xt)Stdt + a{Xt)StdWt, 5*0 > 0 (5) 

where /r : A —>■ R is the drift coefficient and cr : A —>■ (0,oo) corresponds to the volatility. Let Ft be a 
filtration of F satisfying usual hypothesis and right continuous version of the filtration generated by Xt and 
St- Glearly the solution of the above SDE is an Ft semimartingale with almost sure continuous paths. We 
assume that Xt is observed. 
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An admissible strategy is a dynamic allocation to these assets and is defined as a predictable process tt = 
{'^t = ^t <T} which satisfies conditions, given in (A2) below. 

(A2) (i) is square integrable w.r.t St, 

(ii) A(e2) < 00, 

(iii) 3a > 0 s.t. P{Vt > —a, t G [0, T]) = 1 where Vt = ^tSt + StBt. 

It is shown in [5] that if the market is arbitrage free, the existence of an optimal strategy for hedging an J-t 
measurable claim H, is equivalent to the existence of Fdllmer-Schweizer decomposition of discounted claim 
H* := in the form 

= + ( 6 ) 

where J'o, P), = {L^*} o<t<T is a square integrable martingale starting with zero and 

orthogonal to the martingale part of St, and } satisfies A2 (i). Further appeared in the 

decomposition, constitutes the optimal strategy. Indeed the optimal strategy tt = is given by 

it := if, 

Vt* := Ho+ [ iudS: + Lf, (7) 

£t := Vt*-itS:, 

and BtVf* represents the locally risk minimizing price at time t of the claim H. Hence the Follmer-Schweizer 
decomposition is the key thing to verify for settling the pricing and hedging problems in a given market (in 
particular when it is incomplete). We refer to [H] for more details. 

It can be shown in the similar line of [5] that this market model admits existence of an equivalent martingale 
measure. Hence under admissible strategies the market is arbitrage free. To price a claim H of European 
type in the above incomplete market, we would consider the locally risk minimizing pricing approach by 
Follmer and Schweizer, i.e., decomposition of type © and then show that the strategy, thus obtained is 
admissible. 


3 Existence 


Consider the following initial boundary value problem which is known as B-S-M PDE for each i 

d7]i{t,s) 1 2/.N 


dt 






= r{i)r]i{t,s) 


( 8 ) 


for {t,s) G (0,T) X (0,oo) and r]i{T,s) = {s — K)'^, lims4_o ? 7 i(t, s) = 0 for all t G [0,T]. This has a classical 
solution with at most linear growth (see [H])- We also introduce a log normal probability density function 


a(x-, s, i, v) := 


By a direct calculation one has 


-1/32 

e 2^ 




/d(x,s,i,v) := 


ln(f)-(r(*)-^)u 


c^(i)^/v 


+ rU) ^ ^ = 0. 

dv a(i)^ 2 V 2y/v 


(9) 


Set B := ^(f : V ^ R, continuous | lims 4 ,o s, •, •) = 0, ||(/7|| := sup^, | |< ooj . 
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Lemma 3.1. Consider the following integral equation 


(p{t,s,i,y) 


with \im(p(t, s,i,y) 

s^O 


l-F{T-t + y\i) 
l-F{y\ t) 


rT-t 




fi.y + y\ i) 




1- F{y\ i) 


poo 

'y^Pij{y + v) / (p{t + v,x,j,Q)a{x; s,i,v)dxdv 

Jo 


oViG [ 0 ,T],iGx, ye [o,<] 


(10) 

( 11 ) 


where y is as in ([8]). Then (i) the problem (fTOll - dTlT) has unique solution in B, and (ii) the solution of the 
integral equation is in ip(t,s,i,y) is nonnegative. 


Proof, (i) We first note that a solution of (Unil-dll]) is a fixed point of the operator A and vice versa, where 


A(p{t,s,i,y) ■■= 


l-F{T-t + y\t) 


l-FiylO ■’'■<‘•'> + 7 . ' 


r \-r{i)v f(.y + ^\ ^) ^ 


ip{t + V, X, j, 0 ) 0 ( 2 ;; s, i, v)dxdv, 


Clearly B is a. closed subspace of a Banach space {B, || • ||), where B is the set of all continuous functions 
with at most linear growth in s variable. It is also easy to check that for each Lp € B, Aip : (0, 00 ) 

is continuous and limsj^o s, ■, •) = 0. Now in order to show existence and uniqueness in the prescribed 
class, it is sufficient to show that A is a contraction in B. Because, then A : B ^ B is also a contraction 
and Banach fixed point theorem ensures existence and uniqueness of the fixed point in B. To this end, we 
need to show, for (pi,(p 2 G B, \\Aipi — Aip 2 \ \ < J\\ti — 75211 where J < 1. Indeed 


\\Aipi-ALp 2 \\ = sup 

■D 


= sup 

V 


Aipi — Aip2 


1 + s 

pT-t 


-r{i)v fiy±ji\J)_ 


/O 


1 - F{y I i 


J'^Pijiy + '>^) [ (Pi - ‘P2){t + 0) °^^^f'‘^f^'^^ dxdv 


< sup 

V 


= sup 

V 


r + H Vp,^.(y + ^) Hil + x) sup 

Jo 1 - -p^(y I *) ^ Jo V 

Jo ^-F[y\i) I + s 

, / X poo, , s / . N , 

where, a(s) = Jq {1 + x)a(x-, s,i,v)dx = 1 + e V / = 1 + se 

Thus, \\Aipi - Aq} 2 \\ < J\\q}i - 7 ^ 2 II where. 


751 - 752 


I + a; 


1 + s 
a(x] s, i, v) 


I + s 


dxdv 


J = sup 

D 


pT-t 


< sup 


I 


2,"Vi-i"(yl*) 


-r{i)v f{y + v\ i) 1 + 

I - F{y I z) I + s 

fT-t 

f{y + v\i)dv 


dv 


sup 

25 


F{y + T-t\i) - F{y\i) 


< 


1 - F{y\i) 
1-F{y\i) 


1 - F{y\i) 


= 1 


using r{i) > 0 and (Al). 

(ii) Using (AI) and smoothness of rji for each z, the first term on the right hand side is in C^'^’^(T>). 
Under assumption (AI) the fact, the second term is continuous differentiable in y and twice continuously 


5 































differentiable in s, follows immediately. The continuous differentiability in t follows from the fact that the 
term (p{t + v, x,j, 0) is multiplied by ^^((0, oo)) functions in v and then integrated over v € (0, T — <). Hence 
(p{t,s,i,y) is in 

(iii) We note that the problem (|S]) has a closed form solution and that is nonnegative. Hence nonnegativity 
of (fi follows from nonnegativity of all the coefficients of Volterra equation (nni). □ 

Proposition 3.2. The unique solution of (I10D - (I11I) also solves the initial boundary value problem ((¥|) -([ ^ . 


Proof. Let ip be the solutions of (ITTO - (fTT]) . Thus using (fTOj) . ip(T,s,i,y) = T]i{T,s) = (s — K)'^, i.e., the 
condition ([2]) holds. From Lemma lOl fiil. ip is in Hence we can perform the partial differentiations 

w.r.t. t and y on the both sides of (US. We obtain 


dt 


, . f{T-t + y\i) l-F{T-t + y\ i) dp,{t,s) _ fiv + T - t \ i) 

1-F(y|*) {l-F{y\^)) dt l-F{y\^) 


'^Pij{y + T - t) / ip{T,x,j,0)a{x-,s,i,T - t)dx + / e" 




l-F{y\ t) 


^OQ 

^PijiyF'^) / -^(t + v,xJ,0)a{x;s,i,v)dxdv 


( 12 ) 


by differentiating w.r.t. t under the sign of integral. Now, after taking partial derivative w.r.t. y on both 
sides of (nni), we get 


d , . . f(T — t + y \ i) , . 1 — FiT — t + y\i) fiyV) 


dy 


+E 


l-F{y\ i) 

-rWt, F]'{y + v) r 

3 1 - F(?/ I i) io 


ip(t + u, x, j, 0)a(a;; s, i, v)dxdv. 


where q'{u) := -^q{u) and qij{u) = f{u\i)pij(u). For further simplification, we would simplify the last 
additive term in the right side using integration by parts w.r.t. v where q' is treated as second function. We 
would also use the following identities §^ip{t + v,x,j,0) = ^ip{t + v,x,j,0) and ^ = —a to 

obtain 


dy 


Tit,s,i,y) = - 




1- F{y\ i) 


{l-Fiy\t)r 


{p{t,s,i,y) - - 


l-F{T-t + y\ t) 


- F{y I i) 




l-F{y\ t) 

^-r(i){T-t) fjF -tFy\i) 
l-F{y\i) 

fiy I i) 


'^P^J{y + T-t) I ip{T,x,j,0)a{x-,s,i,T -t)dx- Y^-——^'^p,j{y)ip{t,sJ,0) 


- e a(a:;s,f,u)| -r(i)^p„(?/ + u)(^(t + u,a:,j,0) 
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By adding equations (IT^ and (fT^ . we get 

d d 

— s, i, y) + s, i, y) 


l-F{T-t + y I i) dy,{t,s) f{y\i) 

1 — F{y I *) dt 1 — F{y \ i 

f{y + v\i) 

l-Fiy\i) 


rT-t 


- <f{t,s,i,y) - '^p,j{y)(p{t,sj,0) \+ e 


—r{i)v 


{ d 13 1 \ 

+ J ‘p{t + v,x,j,0)a{x;s,i,v)yr{i)+l3-^ + —jdvdx. (14) 


Now we differentiate both sides of (fT0)l w.r.t. s once and twice respectively and obtain 


9 I-F{T-t + y \i) drii{t,s) 

= 1-F(,|,) -* 




/.OO ^ 

/ (p{t + v,x,j,0)a{x;s,i,v) _ dxdv, 

Jo cf{i)^/vs 


(15) 


92 


1 - F{T -t + y\i) d^yi{t, s) 
ds^ 


1- F{y\ i) 

J p{t + v,x,j,0)a{x;s,i,v)-^ 

From equations (ITSl) and (fTBl) . we get 
1 2 ..^ 2d'^P 

r(.).^ + r (■>” ^ 

I - F{T -t + y \ i) ( dr],{t,s) 1 


/3 


1#* 

1 


a‘^{i)v a{i)^/v a‘^{i)v 


dxdv. 


l-F{y\ i) 




a, +r">* a,^ 


rT-t 


-r(i)v fiy + '>J I i) 
l-F{y\i) 


l-F{T-t + y\i) 
l-F{y\ i) 


dyi{t,s) dr]i{t,s) 1 2 /.N 


5t 


+ r(i)i 


a. +r<’> a.^ 


fjy I 

1 - J^(y I i) 


'^P^j{y){+{t,s,j,0) - (p{t,s,i,y)) +r{i) (^{t,s,i,y) - -— J ^ 


fjy I 0 

1-F(2/ I i) 


'^P^Jiy)ip{t,s,j,0) - (p{t,s,i,y)) +r{i)(f{t,s,i,y). 




(16) 


'^P^J{y + v) [ (pit + v,x,j,0)a{x-,s,i,v) ^ - pjlp - dxdv. (17) 

^ Jo \a[%)y/v 2v 2y/v 2vJ 

Finally, from equations (fTOl) . ([8]), ([9]), (fldl) and (13 we get 

d d 9 19^ 

— (p(t, s, i, y) + s, i, y) + r{i)s — <f{t, s, i, y) + s, i, y) 


□ 


Thus equation ([4]) holds. 

From Lemma [3.II and Proposition 13.21 it follows that (jUlhas a classical solution. We prove uniqueness in the 
following section. 


4 Uniqueness 

Proposition 4.1. Assume (Al)(i)-(iii). Let tp be a classical solution of (jd]). Then (i) + solves the integral 
equation (llOD - lfTTll ; (ii) (s — iF)+ < (p{t, s, i, y) < s. 
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Proof, (i) Let {fl, T, P) be a probability space which holds a standard Brownian motion W and the Poisson 
random measure p independent of W. Let St be the strong solution of the following SDE 

dSt = St{r{Xt)dt + a{Xt)dWt), Sq > 0 

where Xt is the solution of ([3]) defined on {Ct,P,P) and driven by p. Let Pt be the underlying filtration 
generated by {St,Xt) satisfying the usual hypothesis. We observe that the process {{St, Xt,Yt)}t is Markov 
with infinitesimal generator A, where 

A<p(s,i,y) = ^{s,i,y) + r{i)s^{s,i,y) + ^a'^{i)s‘^^^{s,i,y)+ ^ X,j{y){(p{s,j,0) - ip{s,i,y)) 

y ^ * jp 


for every function p which is compactly supported in s and in y. If p is the classical solution of (|3|) 
then by using the Ito’s formula on Nt := e~-^o St, Xt,Yt), we get 


dNt = 


= p-f{r{X,,)du 


din 

~'''{Xt)(p{t, St, Xt, Yt) + -Tx{t, St, Xt, Yt) + Mp(t, St, Xt, Yt) 

at 


dt Y dJVlt 


where Mt is a local martingale. Thus from (|4l) and above expression, Nt is also an Pt local martingale. The 
definition of Nt suggests that there are constants ki and fe such that \Nt\ < ki + fc 2 <St for each t, since p 
has at most linear growth. Again, since the following expression 

St = So exp (r(A„) - icr(A„)^)(iu + J a{Xu)dWi 

holds, one concludes that St is a submartingale with finite expectation. Therefore Doob’s inequality can 
be used to obtain Asup^gjQ |A^s| < oo for each t where E is the expectation w.r.t. P. Thus {Nt}t is a 
martingale. Hence 


ip{t,St,Xt,Yt) = ^ Ele^o-XYduj,^^ | ^ Ele-- K)+ \ St,Xt,Yt]. 

By conditioning at transition times and using the conditional lognormal distribution of St, we get 
ip{t,St,Xt,Yt) 

= E[E[e-^"^^^-^‘^-{ST-K)+ I St,Xt,Yt,T^^t)+i] I St,Xt,Yt] 

= P(T„(,)+i > T I Xt,Yt)E[e-f"’^^^-^^^{ST - K)+ \ St, Xt,Yt,T„^t)+i > T] 


+ [ I St,Xt,Yt,T„^t)+i = t + v] 

Jo 


fit - T„^t) +v\Xt) 
1 - F{Yt I Xt) 


dv 


1 - E{T - T^(t) \Xt) , ~ 

^ _ p(Yt I Xt) ^ i 

'^PXtjiYt + v) 


-r(Xt)v f(^t+v\ Xt) 

1 - F{Yt I Xt) 


l-F{T-t^Yt\ Xt) 




\p2fEu{Xt)\fvx 


-dxdv 


fT-t 


1 - F{Yt I Xt) 
+ ^) 


VXtit,St) + 


-r(Xt)v f(Xt + v\Xt) 

1 - F{Yt I Xt) 


' g Pi^t) 2 **)■") 

p(t + V, X, j, 0)- ■;= -—- dxdv. 

' V^xaiXt)y/^ 


(18) 


Finally by using irreducibility condition (Al), we can replace (St, Xt,Yt) by generic variable {s,i,y) in the 
above relation and thus conclude that p is a solution of (nnn-diii)- Hence the proof. 













(ii) We note that 

- K) I Ft] < - K)+ \ Ft] < i?[e"| Ft]. 

Now using martingale property of e~ -^o r{x^)dug^^ from (fT51) and above we get St — K < ip{t, St,Xt, Yt) < St. 
Again, from (1181) . we know that ip is an expectation of a nonnegative quantity, hence nonnegative. Thus (ii) 
holds. □ 

Theorem 4.2. The initial-boundary value problem ^has a unique classical solution in the class of functions 
with at most linear growth. 

Proof. Existence follows from Lemma l3. II and Proposition [3^ For uniqueness, first assume that pi and ip 2 
are two classical solutions of (|3])in the prescribed class. Then using Proposition 14. 11 we know that both also 
solve (nnD-iiii). But from Lemma 13.11 there is only one such in the prescribed class. Hence = <P 2 . D 

Remark 4.1. The above theorem can also be proved in a different manner which heavily depends on the mild 
solution techniques Hg/ and Proposition 3.1.2 oim- Such an alternative approach is taken in m to establish 
well-posedness of a special case of O- The reason for adopting the present approach is that, it enables us 
to establish the equivalence between the PDE and an IE in one go. This in turn suggests an alternative 
expression of partial derivative of the solution. In the next section the importance of such representation is 
explained. 


5 Pricing and optimal hedging 


In this section we consider European call option on the stock dynamics as given in Section 2. The terminal 
claim, associated to that option is {St — K)^. We show that its locally risk minimizing price at time t{< T) 
can be written in terms of the solution of (l4|). We further show that the corresponding optimal hedging has 
an integral representation in terms of the price function. 

Theorem 5.1. Let ip be the unique classical solution of in the class of functions with at most linear 
growth. 


1. Let {f,,e) be given by 
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ds 


Then (f, e) is the optimal admissible strategy. 

2. p{t, St, Xt, Yt) is the locally risk minimizing price of {St — K)'^ at time 0 < t < T. 


(19) 


Proof. Under the market model, the mean variance tradeoff (MVT) process Kt (as defined in Pham et al 
m) takes the following form 



/ p{Xs)-r{X,) 

V 


ds. 


Hence Kt is bounded and continuous on [0,T]. We also know that St has almost sure continuous paths. 
Since, H* G L^(H, F, P) for H = {St — K)~^ we apply corollary 5 and Lemma 6 of [H] to conclude that H* 
admits Fdllmer-Schweizer decomposition m with an integrand satisfying A2(i) and L^ being square 
integrable. Therefore, to prove the theorem it is sufficient to show that 


(a) p{t,St,Xt,Yt) = StBt +ftSt for alH < T; 

(b) P{p{t, St,Xt,Yt)>Oyt<T) = l 

(c) ^(p(t, St,Xt-,Yt-) = Ro + Jo ^tdSf + Lt for all t < T; 
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(d) there exists Fq measurable and Ft measurable Lt such that Lt := E[Lt \ Ft] is orthogonal to 
f* a(Xt)SfdWt i.e., the martingale part of and H* = Hq + ^tdS^ + Lt] 


where is the unique classical solution of (|3]) in the prescribed class and (C,e) is as in (ITOl) . 

From the definition of St in m> (a) follows. In Lemma o and Proposition 14.11 it is shown that (p 
is a non-negative function. Hence (b) holds. Next we show the condition (c) and (d). We apply Ito’s 
formula to S'*, X*, Hj) under the measure P and use ([S]), (|3]) and ([3]) to obtain after suitable 

rearrangement of terms, for all t < T 


= <piO,So,Xo,Yo)+ [ 

Jo 


d(p{u^ Sui 7 ^u— ) 

Fs 


^0 


I ; ^u— ; '2') ; ^U— 9 (^ti— ; ^U— 1 ■^)) 

Jr 

-(p{u, Su, Xu-,Yu-)]p{du, dz) 


where the last integral is w.r.t. the compensator of p. We set 

Lt := j [ [^{u,Su^X^_+h{Xu-,Yu-,z),Y^_-g[X^_,Y^_,z)) 

Jo Jr 

-p{u, Su, Xu-,Yu-)]p{du, dz). 


We note that Proposition id.ll fiil implies that the integrand in above expression is bounded by K. Therefore, 
Lt being an integral w.r.t. a compensated Poisson random measure, is a martingale. Again the independence 
of Wt and p implies the orthogonality of Lt to the martingale part of S^. Thus, we obtain the following 
Follmer-Schweizer decomposition by letting t 'f T., 

bP{St-K)+= p{0,So,Xo,Yo)+ [ ^tdS^+LT. (20) 

Jo 

Thus (c) and (d) hold. □ 

Theorem 5.2. Let p be the unique solution of (|T]). Set 


'ip{t,s,i,y) := 


1- F{T -t + y I i) dr],{t,s) 




1 — F{y I i) ds 

ii((in(f)-(Kb-4^)-)^)^ (ln(f) - (r(i) - 


p{t + v,x,j,0)- 


'/^xsa{i)^/v 


^(^)2 


dxdv 


( 21 ) 


where (t, s, i, y) S V. Then -ip)t, s, i, y) = ■§^p{t, s, i, y). 


Proof. We need to show that if (as in (EU) is equal to Indeed, one obtains the RHS of (EU by 
differentiating the right side of (flUl) with respect to s. Hence the proof. □ 

Remark 5.1. It is well known that in a numerical differentiation, an isolated perturbation gets enhanced 
whereas in a numerical integration the same gets reduced. In ES), the function if, a partial derivative of 
p, is given by an integration involving p. Thus the above theorem essentially provides a robust way to find 
-^pft, s,i,y) using a step by step quadrature method. Further we show that ^ constitutes the locally risk 
minimizing hedging of a contingent claim whose price is given by p. We would like to mention that the 
price of other types of options such as barrier options, composite options, basket options etc can also be 
represented similarly by imposing appropriate terminal and boundary conditions. Thus hedging strategy of 
those would also have similar integral representations. Needless to mention that this finding holds true for 
Markov modulated market model, a special case of the present model. 


10 








6 Conclusion 


To our knowledge, the question of robust computation of hedging for regime switching market is not addressed 
in the literature. In this paper an integral representation of optimal hedging is obtained which leads to a 
robust scheme of computation of hedging. Besides, there is another important contribution. The integral 
equation established here, exhibits how option price in regime switching case depends on that in Black- 
Scholes counterpart. In equation (1^ . we indeed express the optimal hedging using Black-Scholes delta 
hedging and some additional terms. It is also important to note that Ay are the only additional(functional) 
parameter from the Black-Scholes case. Using some nonparametric estimation procedure, as in m , it is 
possible to have a consistent estimate of Ay functions, if the regimes are observed. In a recent paper m 
the convergence of corresponding approximate option price functions is investigated. 
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